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Abstract. We compute the one loop corrections to the pole mass of the gluon in the MS scheme 
in the Landau gauge in both the Curci- Ferrari model and the local composite operator formalism 
with Nf flavours of massless quarks. For the latter we determine an estimate for the gluon mass 
using the effective potential of a local dimension two composite operator and find, for example, 
m gluon = ^'lOAjyfg in Yang-Mills theory. 
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The issue of whether the gluon obtains a dynamically generated mass has been a popular 
topic of investigation in recent years. Following the work of ^ HJ [31 [I] who observed that the 
perturbative vacuum of QCD is unstable, one of the main activities has been on the numerical 
evaluation of the vacuum expectation value of the square of the gauge potential, (^A^). Various 
methods have been used to achieve this ranging from combinations of lattice computations with 
the operator product expansion and instanton considerations, [5J H3 13 [SI H3 EH EEH , to a more 
theoretical approach of the local composite operator formalism of, for instance, |12[ I13[ IT U 115] . 
Moreover, there is evidence from phenomenology that the existence of a gluon mass in the 
range of 500-800 MeV may provide a more accurate explanation of various experimental data. 
Indeed, a valuable summary table of current gluon mass estimates has been given in the article 
by Field, ^5]. Whilst the operator \A 2 l suffers from the immediate objection of not being a 
gauge invariant entity, it has been shown how to relate it to a dimension two gauge invariant 
physical operator, which is the minimization of A 2 ^ over all gauge configurations |121 1171 ITS] .This 
operator, albeit non-local, reduces to a local operator in the Landau gauge and it is solely in 
this gauge that, for example the lattice results of [31 H3 ED have been determined. Indeed the 
local composite operator (LCO) formalism of |12| IT3j was originally developed in the Landau 
gauge but recently an estimate of {\A 2 ^ has been determined in arbitrary linear covariant gauge, 
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Whilst there is much activity in trying to ascertain the existence of a dynamical gluon mass, 
there appears to be less effort into standardizing mass estimates. For instance, in the quark 
sector of QCD the estimates of the various quark masses by methods such as sum rules, lattice 
regularization and the operator product expansion are all expressed as the MS running mass 
at the scale of 2GeV. Although clearly measurements are not always made at this scale. To 
connect the mass estimates one requires an as accurate as possible evaluation in perturbation 
theory of the quark mass anomalous dimension in the MS scheme. This is currently available 
at four loops, ^3 HOI 021 • Moreover, the relation between the quark pole masses and 
the running mass is known at three loops, [21]. For the same problem for a gluon mass the 
analogous quantities are not available to as high an order. For instance, the running of the naive 
gluon mass operator, \A 2 l — acc, in the non-linear Curci- Ferrari gauge, |23J, is known at three 
loops, |2E1- I n the Landau gauge, it transpires that it is not an independent renormalization 
being the sum of the gluon and ghost anomalous dimensions, which is a result that derives 
from a Slavnov- Taylor identity, [221 • This has recently been exploited to obtain the four loop 
running in the Landau gauge for the SU(3) colour group, [2S]. However, the relation between 
the pole mass of the gluon and the running gluon mass is not yet available for QCD at one 
loop. Therefore, it is the aim of this article to provide such a relation for QCD which will build 
on the Yang-Mills expression recently given in ^29, for the Curci-Ferrari gauge. Moreover, since 
the LCO formalism has provided estimates for a dynamically generated gluon mass which are 
comparable with other methods we will also determine the relation for that approach as well. 
This will provide a clean estimate for a gluon mass, since in |12j the effective potential for the 
operator \A 2 L was developed at two loops in the Landau gauge. However, there the estimate for 
a dynamical gluon mass was based on determining the value of an effective gluon mass which 
was by definition a classical mass. It seems to us that a more appropriate quantity to estimate 
through the effective potential approach would be a one loop quantity derived from the gluon 
two-point function such as the pole mass. This is the second aim of the article. 

We begin by defining our notation. We recall that the QCD Lagrangian in a linear covariant 
gauge is 

L QCD = - ^G%G aflu - 7^(d M ^) 2 " c^IV" + i^W 1 (1) 
where a is the gauge fixing parameter, = d^Af, — d u A^ — gf abc A b A c u and i\) tl is the quark 
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field. The indices take the following ranges 1 < a < N A , 1 < I < -/Vp and 1 < i < Nf where iVp 
and iV^ are the dimensions of the fundamental and adjoint representations respectively, Nf is 
the number of quark flavours and f abc are the colour group structure constants. The covariant 



derivatives which determine &t v are 



D^c a 



d^c a 



gf 



abc Ab c 



d^ u + igT a Ay 



>a Aa n /,il 



(2) 



In [T2|, the LCO formalism was derived which involves an additional scalar field a which is 
related to the dimension two composite operator \A 2 . The relevant Lagrangian is 

LCO _ ^QCD _ a _|_ \ aA a A atl — - ^A a A af 



2g 2 C(g) 



——aA a u A a ^ 

2g((g) 



8((g) 



(3) 



where there is an extra contribution to the quartic gluon interaction and L^CO contains the 
usual covariant gauge fixing terms though we will only consider the Landau gauge case, a = 0. 
The quantity £(g) is a function of the coupling constant which has been computed to 0(g 2 ) in 
the Landau gauge in ^21 EC3 El and is such that it ensures the generating functional underlying 
the formalism satisfies a homogeneous renormalization group equation, ^2]- F° r this article we 
note that the relevant terms are 

{\3C A -8T F N f ) 



g 2 C(g) 



9N A 



+ 



(26854MC\T F N f - 1391845Cj - 2l3A<d8C 2 A C F T F N f - 1901760Clr|A^ 2 
+ 221184C ACpTpNf 2 + 58A192C A N^T F - 55296C F TfN? 



65536T N - 

F f J 5184tt 2 N a (35C a - 16T F N f )(l9C A - 8T F N f ) 



(4) 



In j!2l 1131 114j the a field develops a non-zero vacuum expectation value when one computes the 
one loop effective potential of a which is 



V{a) 



9N A 



Aid 



12 



+ 



64 

+ c\ 



,< I 351 ^ x x 351 ^ x x 249 x 27 X 
' ^ _C/ ^ A i A 2 + ^-L/irAiA 3 - y^A 2 + — A 3 



— A1A2 + — A1A3 
lb 62 



+ 



CpAo H CpAa 

128 32 32 



g 2 N A a>* 



+ 0(5 4 ) 



where we have set 
X 1 = [13C A - STpNf}- 1 , A 2 



a 



[35C A - IQTpNf}- 1 

9N A 
[\3C A -8T F Nf) 



[19C A - 8T F N f }- 1 



a 



(5) 

(6) 
(7) 



and \x is the usual MS renormalization scale, which is introduced to retain a dimensionless 
coupling constant in dimensional regularization. 

Now we consider the relation between the pole mass and the running gluon mass in the 
Curci- Ferrari model, |2H], which includes the BRST invariant mass operator 



L 1 



-m 2 A a u A atl 
2 M 



2— a a 

am c c 



(8) 
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where m is the bare mass. With this term the gluon and ghost propagators in the Landau gauge 
are 



(k 2 



It* 



k 2 



5 ab 



(9) 



respectively. With these it is a straightforward exercise to compute the one loop correction to 
the gluon two-point function. In this respect the one loop snail diagram derived from the quartic 
gluon interaction cannot be neglected in the massive case. The result of our computation for 
the pole mass in the Curci-Ferrari model is 



1 + 




128 



C A 



m 2 (//) (10) 



where m = m(fj,)Z m is the bare mass, m(/x) is the running mass and \i is the renormalization 
mass scale. We have renormalized with the usual one loop MS renormalization constants. As a 
check on the expression, we note that it reduces to the same relation given in [22j when Nf = 0. 
Moreover, we have verified the expression of [29; for arbitrary a prior to specifying the Landau 
gauge which provided a non-trivial check on the symbolic manipulation programmes we used 
for this article. 

We have repeated the above computation for the LCO Lagrangian where the bare mass is 
now defined to be a/[g((g)], |T2jEilE], which at leading order is 



(13Ca - 8T F N f ) 
9N A 



9° 



(ii) 



With the additional interactions the expression for the LCO pole mass is of a similar form 



m lco 



1 + 




64 



In 



m 2 (fj,) 



It f n,) jL + 0(g* 



ll7T^ N 

128 , 

m 2 (//) 



C A 



(12) 



for massless quarks in the Landau gauge. Equipped with this result we can now extend the 
method of |121 1131 114j for estimating a gluon mass. In |12l I14j the minimum of the effective 
potential (J5J) was determined by solving dV j£' = 0. Since the factors multiplying the classi- 
cal effective mass are coupling constant dependent, this is equivalent to extremizing V eS (m 2 ). 
However, it seems to us that an alternative approach is to compute instead the extremum of 
yeS^ m 2 ^ w here one inverts (|12|) to obtain m(/z) as a function of m^QQ and then substitutes 
this into ©. Thus we find 

29 
128 

-Cf\\\2 H — tttCf^i^s — ^~r^i 



V 



off 



( m LCo) 



9 



;Ai + 



^c F x 2 + ^C F X 3 



C A 



16 



64 



99 



C\ 
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Ai Ao H Ai A3 

16 32 
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27 
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C A Xi In 



m 



LCO 



7T^ 



(UC A - 8T F N f Y 
81N A 



g C {g)m 



LCO 



(13) 
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Repeating the process to find a minimum necessitates solving 
= 



9 ^ 

2 Al + 



+ C A 



— —C F \\\2 + -r^C F \\\^ - rrrAi 
8 lb 128 

249 . 27 . 99 \ 

- A 2 H A 3 H 7rV3Ai 

128 64 128 V 

81 81 
A1A2 + ttt;AiA3 



16 



32 




27 / 

which corresponds to the condition 



m 



LCO 



7T^ 



(13CU - 8?>iV^ 
81iV A 



■5 2 C 2 (ff)"iLCO 



(m 2 



LCCV 



dm lco 



. 



(14) 



(15) 



We have not in fact substituted for the explicit expression for Q(g) since this function factorizes 
off the expression for the location of the minimum. If we were to include that part of the 
series which was already known it would introduce an unnecessary truncation error into our 
final estimates for the pole mass. At this point to solve for the mass a scale needs to be chosen 
for [i. In |12l 1131 114j . the choice of scale was such that it removed the logarithm terms. For 
this potential we will take a more general approach and instead set "t-lco = s ^ 2 where s is an 
arbitrary parameter. This means we have determined an equation for the value of the coupling 
constant function of s. In other words 



y 



36C A (16T F N f - 35C A ) [(3465^3 + 4620 ln(s) - 25690) C\ - 864C F T F N f 

+ (19240 - 1584vr\/3 - 3792 ln(s)) C A T F N f 



+ (768 ln(s) - 3328)7$^ 



-1 



(16) 



where y = C A g 2 /(167r 2 ). Through the definition of the running coupling constant we have the 
one loop relation 



where 



g 2 M 
16vr 2 



00 



A) In 



A 2 — 

L MS J 



11 



-c A 



-T F Nf . 



(17) 



(18) 



Hence, we can relate the coupling constant to the scale n and A^jg and deduce a value for n-LCO' 
We find 



m LCO 



A« exp 



((346571-^/3 - 25690) C\ - 864C F T F N f 
+ (19240 - 158471-^) C A T F Nf - 3328T£iVy) 

24 {11C A ~ 4T F N f ) (35C A - 16T F N f ) 



(19) 



which is the main result of this article. It turns out that this is independent of the intermediate 
parameter s. In other words, no matter what scale \i one chooses, one will always obtain the 
same value for the solution to (j!4j) at one loop. 
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N f 


™SU(2)/^ 


™SU(3)/*m 





2.10 


2.10 


2 


1.54 


1.74 


3 


1.24 


1.55 



Table 1. One loop estimates of the gluon effective mass for SU(2) and SU(3). 

We have given the explicit values of the pole mass estimates from (fT9|) for SU (2) and SU (3) 
in Table 1 . Compared with the classical effective gluon mass estimates of |121 1131 114j the Yang- 
Mills estimates have increased by about 5% for SU(3). However, for Nf ^ there is a significant 
decrease. Although this is disappointing it is important to recognise that since they have been 
derived in a scale independent and therefore renormalization group invariant manner, they may 
be closer to the true result, though the inclusion of quark mass may alter these estimates. 

We conclude with several remarks. First, we have constructed a one loop renormalization 
group invariant pole mass for the gluon using the LCO effective potential of |121ll3llT4*| . However, 
it would be interesting to see whether this feature persists at the next order. This only requires 
an extension of the present one loop result since the two loop LCO effective potential is available. 
Although we have ignored quark mass effects it seems that if one could include quark condensates 
in the LCO formalism in addition to that for ^A a ^A'j 1 then it might be possible to ascertain 
the extent to which condensates could be responsible for the quark and gluon masses. If the 
renormalization scale invariance persists even at one loop for this scenario then one would 
not have to worry about solving a multi-scale type renormalization group equation. Our final 
comment concerns the situation where a gluon mass is dynamically generated through, say, the 
LCO formalism. If this is the case then one would have to include additional contributions due 
to a gluon mass to the existing quark pole mass multi-loop estimates. 
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